Decay of a charged scalar and Dirac fields in the Kerr-Newman-de Sitter background 
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We find the quasinormal modes of the charged scalar and Dirac fields in the background of the 
rotating charged black holes, described by the Kerr-Newman-de Sitter solution. The dependence of 
the quasinormal spectrum upon the black hole parameters mass M, angular momentum a, charge 
Q, as well as on values of the A-term and a field charge q is investigated. Special attention is given 
to the near extremal limit of the black hole charge. In particular, we find that for both scalar and 
Dirac fields, charged perturbations decay quicker for g > and slower for q < Q for values of black 
holes charge Q less than than some threshold value, which is close to the extremal value of charge 
and depend on parameters of the black holes. 

PACS numbers: 04.30.Nk,04.50.-|-h 



I. INTRODUCTION 

Quasinormal spectrum of black holes has been exten- 
sively investigated during recent years for a great variety 
of black hole backgrounds and fields, because it is an im- 
portant^characteristic for observation of the gravitational 
waves [li], stability analysis AdS/CFT calculations of 
temperature Green functions Special attention has 
been paid to perturbations of a scalar field [J, as a sim- 
plest model, when the influence of the spin of the field is 
neglected. When one considers the charged black hole, 
the scalar electrodynamic can successfully model the in- 
teraction of the charged field with the electromagnetic 
background of the black hole. Therefore the calcula- 
tion of the quasinormal modes of charged fields, initiated 
in for charged scalar field in the Reissner-Nordstom 
and dilaton backgrounds, was continued in the further 
research j6|, lS, i^, [!]■ In particular, in Q the quasi- 
normal modes of the massive charged scalar field were 
found with the WKB accuracy. In [| and R it was 
shown that the quasinormal modes, corresponding to the 
charged scalar field, decay quicker than those of the neu- 
tral field unless the black hole charge is larger than some 
large near extremal value. This is opposite to the behav- 
ior at asymptotically late times, characterized by the so- 
called "tail" decay, when the charged scalar field decays 
slower, and therefore, dominates at asymptotically late 
times Yet, the quasinormal frequencies calculated 

for Reissner-Nordstrom black holes in [5, ^, JJ with the 
help of the WKB method need better accuracy and can- 
not be trusted near the extremal limit, especially for the 
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scalar case, because the effective potential is frequency 
dependent in this case, and WKB equation for the QN 
frequency must be solved together with a frequency de- 
pendent equation determining the position of the maxi- 
mum of the effective potential. 

On the other hand, one has much richer physical sit- 
uation, when one takes into consideration all the rele- 
vant parameters, such as black hole angular momentum 
and the cosmological term, i.e. when one starts from the 
Kerr-Newman-de Sitter (KNdS) black hole as a gravita- 
tional background. In this paper we achieve both the 
above mentioned aims: to find QN modes with very high 
accuracy, by using the convergent Frobenius expansion, 
and to consider the most general relevant black hole solu- 
tion of general relativity, KNdS solution. The latter gives 
us dependence of the QN spectrum on a great number of 
parameters: charge of the black hole Q, charge of the 
field q, normalized angular momentum of the black hole 
a, the cosmological term A. 

The paper is organized as follows: Sec. II gives some 
basic formulas on KNdS metric and on radial wave equa- 
tion for charged massless scalar and Dirac fields, and also 
discusses the system of units we used for showing the 
QNMs. Sec. Ill reviews the obtained numerical results. 
In the Conclusion we summarize the obtained results. 



II. BASIC FORMULAS 

In the Boyer-Lindquist coordinates the Kerr-Newman- 
de Sitter metric has the form [Till 



Ae sin^ 9 
(1 + a)2p2 



[adt - (r^ +a'^)dip]'' + 



(1 



(dt — a sin^ OdipY , 



(1) 



2 



where 



= (r2 + a2) (1 - arVa^) - 2Mr + Q2 ^ = AaVS 
Ae = 1 + a 008^6', 



p2 — 7-2 _j_ q2 j,pg2 I 



The electromagnetic background of the black hole is given by the four- vector potential 

Qr 



Aadx'' = — ^ [dt - asin^ dd^p). 



(2) 



(3) 



(1 + a)2p2 

The charged scalar and Dirac fields in curved space-time are described by the following equations of motion : 

^fX^^" - iqA^.g^'" [2^^!^) - - iqA^,,^g^'' ^'^^^ = 0, (charged scalar) (4) 



y"e^(5^ + r^, + = 0, (charged Dirac) 



(5) 



where q is the charge of particles and Aa is the electromagnetic potential of the background. 

Some manipulation with scalar and Dirac equations allow in the some separable form (see [ll| and references 
therein). The existence of the Killing vectors dt, d^, implies the exponential harmonics of the following form ~ e~'"*, 
~ e""*^). After the separation of the angular, radial and time variables 

<^'-'^t,r,e,(t)) cx e-*"*e™*S'(6l)i?(r), 

one can obtain the equation for the angular part 

(1 + afw? 



d 9, ,^ os d . ,^ , (1 + a)^ 2 



dx dx 



1 + a 
1 + ax'^ 



2s{am — (1 -|- Oi)(0x 



(1 + a) 



(1 + a){s^ + 2smx) 



(1 + aa;2)(l - x^) 



-C^ - 2m(l + a)C + 



>Ss{x)=Q 



(6) 



where s is the spin of the field (s = 0, 1/2 for the scalar and Dirac field respectively), x = cos0, ^ = aw, and A 
is the separation constant (for the non-rotating case A = ^{l -f 1) — s{s — 1), where ^ > s is positive (half)integer 
multipole number). The angular equation can be solved numerically with respect to A for each value of lo by using 
the three-term recurrence relation illil. 
The equation for the radial part is [ll| 



1 



A-^— A^+i — + - 
^ dr ^ dr A 



dAr 



2a 



K-^]+4:is{l + a)ujr- — {s + l){2s + iy + 2s{l-a)-2isqQ-X }Rs{r)^0, (7) 



where K = [Lu{r^ + a^) — am]{l + a) — qQr. 



r 



Generally, this equation has five regular singularities: 
the event horizon r+, the internal horizon r_, the cos- 
mological horizon r^, the spatial infinity and one more 
singularity at r'_ . One should note, that r± and r'j_ are 
roots of equation A,. = 0. 

In the limit of A ^ one has, 



r±^ M± y/A'P -a^-Q"^, 



±- 



/a 



and, therefore, the spacial infinity becomes an irregular 
singularity 

By the definition, the QNMs are eigenfrequencies uj 
of (O which satisfy the following boundary conditions 
(b.c): 'ijj represents purely ingoing waves at the event r = 
r+ horizon and purely outgoing waves at the cosmological 



horizon r — r'^ (or the spacial infinity if A = 0). 

Now we introduce the new function, which is regular 
at these two points if the QNM b.c. are satisfied 



y[r) = r 
dB{r 



,2s+l 



s+2iX(r+)/A;,(r+) 



*^('^)i?(r). 



where 



dr 



K 

Ar 



(8) 



The appropriate Frobenius series are 

f r - r+ Y f 1 - pr-/7 

where p = r^/r'^. 



1-p 



(9) 



3 



Substituting ([8]) and ([9]) into one can obtain the 
A'^-term recurrence relation for the coefficients a, 



'min{N — l^i) 



J=0 



3,1 



(lo) a,-j =0, for i > (10) 



where N depends on the black hole parameters. For the 
particular Schwarzschild case iV = 3, but for the general 
case under consideration N is higher than 3. We decrease 
the number of terms in the recurrence relation using the 
Gaussian eliminations: 



c},/(^) 



(a;) 



for J = 0, or i < k, 



After one finds / numerically, he can solve the equation 
with infinite continued fraction (see [l^ for more details) 



„(3) 



„(3) „(3) „(3) (3) 
^2,n'--0,n-l ''2,n-l''0,n-2 



„(3) 



„(3) 



„(3)J3) 

^1,1 



(3) ^(3) ^(3) ^(3) 



X),n+l'^2,n+2 '^0 ,n+2'^2 ,n+3 



Co 



-'l,n+2 



.(3) 



(11) 



Since we can find the separation constant A for each 
particular value of w, pT|) allows to find QNMs with the 
desired precision. This technique of the QN spectrum 
calculation was proposed by Leaver [13]. 

Now we shall discuss the units of measurements and 
ranges of the black hole parameters. In this paper we 
shall measure all the quantities in units of the event hori- 
zon. For this we choose the black hole mass, so that 
A,.(l) = 0. Then one has, 

2M= (^1-|) (l+a') + Q'. 

We parameterize the cosmological constant A by the pa- 
rameter p < 1 



A, (1/p) = 



A 



1 - p(a2 + Q2) 



3 ^ l + p + /92(a2 + Q2)- 



In these units the condition AJ,(1) > gives us the range 
of values of the black hole charge. 



l + 2p 



The positivity of the righthand side of the above equation 
bounds the possible values of the parameter of rotation 




FIG. 1: Charged (e) scalar field fundamental frequency {I = 
0) for the Reissner-Nordstrom black hole as a function of its 
charge (Q). q = —0.3 (purple), q = —0.2 (magenta), q — —0.1 
(red), q — (orange), q = 0.1 (green), q = 0.2 (cyan), q = 0.3 
(blue) . The larger q corresponds to the larger real and larger 
imaginary (for small black hole charge) part of the QNM. 





FIG. 2: Charged (e) scalar field fundamental frequency {I = 
0) for the Kerr-Newman black hole of charge Q — 0.2 (a < 
2v^/5) and Q = 0.8 (a < 0.6) as a function of a. q — —0.3 
e — (orange), q — 0.1 (green), q = 0.2 (cyan), q = 0.3 
(blue) . The larger q corresponds to the larger real and larger 
imaginary (for small black hole charge) part of the QNM. 
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TABLE I: QNMs of the scalar perturbation of the Kerr-Newman black holes {a — 0.4). 



£ 


m 


g = 0.2, g = -0.3 


Q = 0.2, g = 


Q = 0.2, g = 0.3 


Q = 0.8, g = -0.3 


Q = 0.8, g = 


Q = 0.8, g = 0.3 








0. 153622 -0.158460i 


0.193444-0.1647774 


0. 234254-0. 169836i 


0.027107-0.085597i 


0.141475-0.105014* 


0.274019 -0.093648i 


1 
1 
1 


-1 

1 


0.396131-0. 156125i 
0.479121-0.151767i 
0.603295-0.1459514 


0.428770-0.160382* 
0.514846-0. 154723i 
0.642846-0.147968* 


0.462269-0.164345* 
0.551295-0.157438i 
0.682964-0.149826* 


0.236903-0.098374i 
0.285439-0.093609* 
0.364922 -0.075864i 


0.337666-0.109255* 
0.399817-0.098162* 
0.505807-0.068374* 


0.452137-0.115708i 
0.528744-0.097728i 
0.667302-0.058104* 


2 
2 
2 
2 
2 


-2 
-1 

1 
2 


0.656529-0.156419i 
0.726943-0.154996i 
0.815194-0. 152025i 
0.927889-0.147965i 
1. 072198-0. 144707i 


0.688008-0.1591792 
0.760246-0.1572874 
0.850587-0.153847* 
0.965608-0.149382* 
1.112314-0.145847* 


0. 720043-0. 161832i 
0.794064-0.159477i 
0.886446-0.155578i 
1.003731-0.150726i 
1. 152765-0. 146933i 


0.443727- 0.102703i 
0.488765 -0.100428i 
0.546533-0.095406* 
0.625000-0.085165* 
0.742952 -0.064746i 


0.543278-0.109491* 
0.595790-0.105290* 
0.663393-0.097610* 
0.756010-0.083476* 
0.897860-0.058622* 


0.651446 -0.114676i 
0.711645-0. 108424i 
0.789555-0.097946i 
0.897414-0.079871i 
1.064754-0.052886* 



TABLE II: QNMs of the scalar perturbation of the Kerr-Newman black holes [Q = 0.2). 



e 


m 


a = 0.2, g = -0.3 


a = 0.2, g = 


a = 0.2, g = 0.3 


a = 0.8, g = -0.3 


a = 0.8, g = 


a = 0.8, g = 0.3 








0.167766-0.184220* 


0.211875-0. 192171i 


0.257011-0.198680i 


0.105432-0.103623* 


0.132795-0.106936* 


0.160904-0.109273i 


1 
1 
1 


-1 

1 


0.470255-0.176446* 
0.522144-0.175672* 
0.585319-0.175043* 


0.507611-0.180941i 
0.561332-0.179498i 
0.626523-0. 178296i 


0.545869-0.185127* 
0.601338-0.183046* 
0.668457-0.181312* 


0.266598-0.109174* 
0.351193-0.097777* 
0.545520-0.054487* 


0.289229-0. 112303i 
0.377110-0.099006i 
0.578091-0.054264i 


0.312498-0.115211* 
0.403534-0.100059* 
0.610920-0.054055* 


2 
2 
2 
2 
2 


-2 
-1 

1 
2 


0.787475-0.176507* 
0.834779-0.176198i 
0.887964-0.175702* 
0.947908-0.175148* 
1.015500-0.174741* 


0.823858-0.179350i 
0.872296-0.178790i 
0.926678-0.178052i 
0.987872-0.177273i 
1.056747-0.176665i 


0.860824-0.182078* 
0.910368-0.181272* 
0.965915-0.180299* 
1.028326-0.179304* 
1.098451-0.178503i 


0.438203-0.109596* 
0.505143-0.106753* 
0.600747-0.097893* 
0.754274-0.075327* 
1.039415-0.048194* 


0.459854-0. 111675i 
0.528641-0. 108258i 
0.626685 -0.098644i 
0.783768-0.075164i 
1. 072980-0. 048219i 


0.481918-0.113676* 
0.552512-0.109685* 
0.652951-0.099329* 
0.813525-0.074974* 
1.106653-0.048245* 



III. RESULTS 

The quasinormal frequencies of the Kerr-Newman-de 
Sitter black holes depend on a number of parameters: 

1) black hole parameters: mass A/, charge Q, angular 
momentum a, 

2) field parameters: field charge q and spin s, 

3) cosmological constant A, 

4) numbers of modes in the spectrum: the multi-pole 
number £, the azimuthal number m, and the over- 
tone number n. 

Therefore, if one wants to represent quasinormal fre- 
quencies for all values of the above parameters, he has 
to show a great amount of table data. We decided to be 
limited here by representative tables or plots, which will 
show dependence of the quasinormal modes on each of 
the above parameters. Thus, for example on Fig. 1 one 
can see the dependence of the Re{uj) and Im{Lu) on the 
black holes charge Q for fundamental mode £ = n — 
of perturbations of the scalar field. The real part of cu 
monotonically grows with the black hole charge Q and 
the field charge q (note that q can be both positive and 
negative). For positive values of q, Re{uj) attains some 
maximum value as a function of Q, at some large value 
of Q. Imaginary part has more complicated behavior: 
it monotonically decreases as a function of Q until some 
near extremal value of the black hole charge, keeping 
meanwhile monotonic dependence on q. Then the curves 
with different q intersect, that is, the larger q does not 
guarantee larger Imuj. Thus if for not very large Q, the 



charged field decays quicker than the neutral one, for the 
near extremal Q, the charged field decays slower than 
neutral. From Fig. 1 one can see that this happens at 
Q ^ 0-SQextr in the units of the event horizon, while 
in "Af — 1 units" this corresponds to Q ~ 0.995Qextr- 
Even though the WKB technique, developed until higher 
orders fl^. reproduces this intersection shown in Fig 1., 
it could not be easily trusted in this case because of uj 
dependency of the effective potential Therefore, con- 
firmation of the intersection with the help of the conver- 
gent and accurate Frobenius method leaves apart possible 
interpretations of charged quasinormal modes in the con- 
text of universality of the critical collapse . Let us note 
also that the non-monotonic behavior for some curves for 
real and imaginary parts of w, near the extremal values 
of charge, depends on the value qQ and is not new in 
fact. When approaching the extremal limit of values of 
Q closely enough, one has the picture of spiraling of the 
plot of Re{Lu)-Im{uj) 

An important point is the checking of all the known 
particular limits for our calculations. For q — Q = 
we reproduce the quasinormal modes for the scalar and 
Dirac fields for Kerr-de Sitter black holes, while for 
A = 0, a = we obtain the limit of pure Reissner- 
Nordstrom black holes and the results of Q, 0. For 
g = we find the quasinormal modes for neutral fields 
for KNdS black holes. When approaching the limit of ex- 
treme values of the A-term, one can reproduce the exact 
quasinormal modes obtained as a solution of the Poschl- 
Teller equation (ITij . if g = a = 0. 

In Fig. 2 we demonstrate the dependence of Re{uj) and 
Im(uj) of scalar field modes on the angular momentum 
a for a few fixed values of the charges Q and q. There 



5 




FIG. 3: Real and imaginary parts of the fundamental QNMs 
of Kerr-Newman-de Sitter BH s = / = 1/2, m = -1/2, Q = 
0.2, a = 0.8, g = (red) (p < 0.51) s ^ I = m = 0, Q ^ 0.5, 
a — 0.5, q — (blue) (p < 0.94 ) as a function of p. 



larger values of the multi-pole number £, under the same 
value m, imply larger Re{ijj) and smaller Im{ijj). Under 
some fixed I, large azimuthal numbers give larger Re{uj) 
and smaller damping rates. When the rotation is large, 
this decreasing of Im{uj) for larger m is considerable, so 
that high m modes can be many times longer lived than 
those for non-rotating case a — 0. This happens for all 
values of Q and q. Finally, on Fig. 3, one can see the rep- 
resentative dependence of the quasinormal frequencies on 
the values of the A-term. The behavior is qualitatively 
the same as for the ordinary Schwarzschild-de Sitter black 
holes [l^, that is, the increasing A-term suppresses con- 
siderably the Reuj and Imuj, independently on the values 
of other parameters. Let us note that when all the pa- 
rameters a, Q and q are non-vanishing, and the A-term 
approaches near the extremal limit, the Frobenius series 
converges very slowly, so that we could not reach there 
extremal limit very closely. 



IV. CONCLUSION 



one can see that for not large values of qQ, one has the 
monotonic decay of both real and imaginary parts of uj as 
a function of a. In other words, the quicker rotating black 
holes have longer lived modes and smaller real oscillation 
frequencies. This monotonic behavior breaks for larger 
values qQ, as it can be seen from the curve Q — 0.8, 
q = 0.3 on Fig. 2. 

In Tables I, III, the quasinormal frequencies for dif- 
ferent values of multi-polar i and azimuthal m numbers 
are given for the first mode n = 0, for some values of 
the charges Q and q, and for a fixed value of the angular 
momentum. In Tables II, IV, on the contrary, we fix the 
values of the black holes charge Q, and consider different 
values of a and q. From the above table one can see that 



In this work, with the help of an accurate convergent 
Frobenius method, we performed an extensive calcula- 
tion of quasinormal modes of charged scalar and Dirac 
fields for Kerr-Newman-de Sitter black holes and have 
analyzed the dependence of the QN spectrum upon the 
great variety of parameters of the black holes Q, M, a, of 
A-term, and of the field parameters q and s. This gener- 
alizes a number of previous works when only some of the 
parameters were considered non-vanishing. The model 
we considered might be successful, when considering the 
interaction of the charged fields with the electromagnetic 
background of rotating black holes. 
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TABLE III: QNMs of the Dirac field perturbation of the Kerr-Newman black holes {a = 0.4). 



e 


m 


Q = 0.2, q = -0.3 


Q = 0.2, q = 


Q = 0.2, q = 0.3 


Q = 0.8, q = -0.3 


Q = 0.8, g = 


Q = 0.8, q = 0.3 


1/2 

ift 


-^1/2 
1/2 


0.246256-0.1541334 
0.342822-0.141572* 


0.280565-0.160311i 
0.380628 -0.144882i 


0.316133-0.165779i 
0.419348-0.147793* 


0.118776-0.090357* 
0.171964-0.078214* 


0.221025 -0.107204i 
0.297229 -0.073804i 


0.344150-0. 113238i 
0.454538 -0.058993i 


3/2 
3/2 
3/2 
3/2 


-3/2 
-1/2 
1/2 
3/2 


0.515027-0.155523* 
0.588694-0.153027* 
0.687238-0.148338* 
0.820628-0.143438* 


0.546793-0.158973i 
0.622736-0.155715i 
0.723951-0.150318i 
0.860126-0. 144895i 


0.579262-0.162249* 
0.657409-0.158246* 
0. 761209-0. 152174i 
0.900057- 0.146265i 


0.332623-0.100337* 
0.378403-0.097191* 
0.442312-0.088796* 
0.543282-0.067958* 


0.432175-0. 108952i 
0.487664-0. 102514i 
0.565845 -0.08904H 
0.693191-0.060714i 


0.542690-0.114981* 
0.608266-0.105015* 
0.701925-0.086213* 
0.858738-0.053488* 



TABLE IV: QNMs of the Dirac field perturbation of the Kerr-Newman black holes (Q = 0.2). 



£ 


m 


a = 0.2, q = -0.3 


a = 0.2, g = 


a = 0.2, g = 0.3 


a = 0.9, g = -0.3 


a = 0.9, g = 


a = 0.9, g = 0.3 


1/2 
1/2 


-1/2 
1/2 


0.287309-0.174116* 
0.340964-0.170124* 


0.325875-0.180887i 
0.381353-0.175322i 


0.365774-0.186924* 
0.422875-0.179958* 
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